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ABSTRACT 



We analyze 3.5 years of public Fermi/~LAT data around the position of the supernova remnant HB21, where four 
point-like sources from the 2 nd Fermi /LAT catalog are located. We determine that the gamma-ray source is produced 
by a single extended source. We model the observed morphology as a uniform circle. The spectral energy distribution is 
best described by a curved power law, with a maximum at 413 ±11 MeV. We divide the circle into three regions defined 
by previously identified shocked molecular clouds, and find that one of these regions has a softer spectrum. The > 3 GeV 
gamma-ray emission of the soft spectrum region is bow-shaped and coincident with the supernova remnant shell seen 
at radio wavelengths. We suggest that the gamma-ray emission from HB 21 can be understood as a combination of 
emission from shocked/illuminated molecular clouds, one of them coincident with the supernova remnant shell itself. 
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1. Introduction 



Gamma-ray emission from SNRs is key for understand- 
ing the processes that lead to particle acceleration in the 
Galaxy and, in particular, to the production of cosmic rays. 
Since protons cannot be traced back to distant sources, pro- 
ton acceleration in SNRs may be probed by observing their 
effects on dense molecul ar clouds located in the vicinity of 
the accelerating region (|Gabici et al.ll2009t) . In particular, 
the observation of gamma rays from neutral pion decay is 
considered to be the smoking gun of cosmic-ray produc- 
tion, since they trace the collisions of accelerated protons 
with nucleons from the ambient medium. Examples of such 
interaction are firmly established in cases like that of the 
SNR W44 (lUchivama et alj|20ll . where GeV gamma rays 
emanate from regions clearly offset from the SNR shell. 

HB21 (G89.0+4.7) is a 19000 year olcQ 
(jLeahv fc Aschenbachl[P996T) mixed-morp hology supernova 
remn ant (SNR) at a distance of 0.8 kpc (jTatematsu et al.l 
Il990t) . As seen in radio continuu m images, the SNR 
displays an elliptical shell of 2° x 1.5° (jCondon et al.lll994l) 
(mean diameter of ~ 25 pc), slightly tilted in the NW-SE 
direction. Only weak, center-filling X-ray emission of 
therm al origin is associated with HB 21 (|Lazendic fc Stand 
12006ft . 

The interaction of HB 21 with the surrounding inter- 
stellar medium (ISM) has been intensiv ely studied. Give n 
the absence of OH masers around HB 21 (|Frail et al.lll996|) . 



Send offprint requests to: I. Reichardt (ignasi@ifae.cat) 

1 We quote here the commonly accepted age for this object, 
but we note that there are indic ations that HB 21 could have 
an age of the order of 5000 years (jLazendic fc Slanell2006r i 



the evidence of interaction between the blast wave and 
the ISM is established by means of local dynamic effects 
that broaden emission lines. In the following paragraph we 
briefly describe the clouds that present such b road emission 
lines. We refer to Fig. 1 in lBvun et al.l ((2006) for a detailed 
view of the interacting clouds, or Fig. [I] in this work for a 
schematic view. 

Evidence of shocked molecular gas was found by 
iKoo et all (|200ll ) in the northern and in the southern parts 
of the shell. The northern cloud (cloud N hereafter) con- 
sists of several small, bright clumps plus a diffuse com- 
ponent extending to the E. The southern cloud (cloud S) 
presents a complex filamentary structure, with a velocity 
spread of up to 40 km s _1 for some particular clumps, and 
it i s coincident w i th a m ass of shocked atomic gas detected 
bv lKoo fc Heilesl (|1991h . There is also a bow-s haped cloud 
at th e NW rim of the radio shell (cloud NW) (|Bvun et al.l 
120061 ). and the central thermal X-ray bright area is occupied 
by small eva porating clouds. Moreov er the so-called clouds 
A, B, and C (jTatematsu et al.lll990T) are aligned N-S in the 
approximately straight E rim of the SNR. These clouds may 
be regarded as overdensit ies of the giant molecular c loud of 
the Cyg OB7 association (jHuang &: Thaddeusll 1986th which 
provides the distance estimate for HB 21. Clouds A, B, and 
C arc located where the eastward blast wave apparently 
collides with the so-called wall. The wall consists of a sharp 
edge of otherwise smoothly distributed atomic gas, which 
extends beyond the SNR boundary, both N and S, therefore 
suggesting that it is a pre-existing structure that affects 
the evolution of the SNR and not the other way around. 
Probably, the wall is the border of the cavity resulting from 
a former HII region around the HB21 progenitor, which 
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might be a former me mber of the Cyg 0B7 association 
(jTatematsu et al.lll990f ). There is also the possibility that 
the coincidence of the A, B, C clouds and the wall with 
the SNR shell is a projection effect, where the clouds are 
indeed in t he vicinity of th e SNR, but they lie in front or 
behind it (|Koo et alJkOQlD . iBvun et all (|2006D suggested 
that the SNR could be as far as 1.7kpc, in which case the 
whole Cyg OB7 complex would be in the foreground. 

According to lNolan et al.l (|2012l) three point-like sources 
(2FGL J2041.5+5003, 2FGL J2043.3+5105 and 2FGL 
J2046. 0+4954) in the second Fermi /LAT source catalog 
(2FGL catalog hereafter) are coincident with the extended 
radio emission of HB 21. In this article we report a detailed 
analysis of the public Fermi /LAT data that will lead to a 
deeper understanding of the gamma-ray emission from this 
object. 

2. Data analysis 

We analyzed Fermi /LAT data with the LAT analysis soft- 
ware, the ScienceTools version v9r27pfl We analyzed Pass- 
7 data corresponding to the period between August 4* /l 
2008 (start of science operations) and February 2, 2012. 
Since HB 21 is almost 5° off the galactic plane, it allows 
for a notably background-reduced analysis with respect to 
the sources found at lower galactic latitudes. To exploit 
this advantage, we defined a region of interest (ROI, i.e. 
the sky region whose LAT photon events are considered) 
as a circle of 10° radius centered on the position (a, S) = 
(20 ft 41 m 05 s , 51°15'58"), which is displaced by 1° toward 
positive galac tic latitudes with respect t o the catalog posi- 
tion of HB 21 (|Ferrand fc Safi-Harbll2012l) . We selected class 
2 events in the energy range between 100 MeV and 100 GeV. 
We applied a set of quality cuts, including the require- 
ment for the spacecraft to be in normal operation mode 
(LAT_CONFIG==l), data to be flagged as good quality 
(DATA_QUAL==1), and a cut on the rocking angle of the 
spacecraft (ABS(ROCK_ANGLE)< 52°). In addition, we 
applied a zenith angle cut of 100° in order to prevent event 
contamination from the Earth limb. Data are binned in sky 
coordinates with the gtbin tool, using square bins of 0.125° 
side. We refer later to the two-dimensional histograms re- 
sulting from gtbin as count maps. 

To study the morphological and spectral properties of 
HB21, we performed a three-dimensional (two spatial di- 
mensions plus the energy) maximum likelihood analysis, 
using the standard gtlike tool. In this method, the likeli- 
hood is computed for different models defined by the posi- 
tion and morphology of the sources producing gamma rays 
in the ROI. For each source, a different spectral shape is 
assumed, and the spectral parameters are left free in the 
likelihood maximization. Our starting point consists of the 
standard galactic and extragalactic diffuse emission models 
provided in the ScienceTools, plus the point-like sources in 
the 2FGL catalog lying up to 15° away of the ROI cen- 
ter, excluding 2FGL J2041. 5+5003, 2FGL J2043.3+5105, 
2FGL J2046.0+49 54 and 2FG L J2051.8+5054. We note 
that according to iNolan et al.l (|2012l ). the source 2FGL 
J2051. 8+5054 is not associated to HB21, but it lies very 
close to the NE edge of the SNR shell, in remarkable co- 
incidence with the above-mentioned cloud A (see Fig. ITa|). 

2 data and software are publicly available at 
http : ll± ermi . gsfc.nasa.gov/ssc 



For this reason, in the following section we consider 2FGL 
J2051. 8+5054 as part of the gamma-ray emission related to 
HB 21. We refer to this model as the null hypothesis, since it 
assumes that the gamma-ray emission from HB21 is faint 
enough to not be distinguished from the background. Next, 
we compare the maximum likelihood obtained with the null 
hypothesis to those of several models that include vari- 
ous morphological descriptions of the GeV emission from 
HB21. First, we include the four point-like sources men- 
tioned above. In addition, we explore the possibility that 
the observed emission is from an extended, resolved source. 
For this, the four point-like sources associated to HB 21 
are replaced by extended source templates. The scenarios 
considered are detailed in sect. 13.11 For each model, the 
goodness of the likelihood fit is estimated by means of a 
test statistic (TS) defined as 

TS = -21og(£ /£), (1) 

where Co and C are the likelihood of the null hypothesis 
and the tested models respectively. Since TS is a likelihood 
ratio of two nested models, it asymptotically follows a \ 2 
distribution with a number of degrees of freedom equal to 
the extra number of par ameters of the test hypothesis with 
respect to the null one jAbdq et al.ll2009al) . Despite some 
caveats (jProtassov et al.ll2002j) . it is normally accepted TS 
= 25 as a detection threshold of a source with two spectral 
parameters (flux normalization and spectral index), which 
corresponds to a statistical significance of 4.6 sigma. In all 
cases, gtlike is run in a two-step procedure: first allowing 
a loose tolerance up to 10% in the fit parameters and using 
the MINUIT method and, second, using the output of MINUIT 
as the initial value for a refitting of the model parameters 
with a tighter requirement of 0.1% tolerance and using the 
NEWMINUIT method. The output of gtlike allows us to gen- 
erate synthetic maps with the expected source shape and 
brightness, given the model best-fit parameters. These syn- 
thetic maps can be subtracted from the counts maps in 
order to visualize the disagreements between the real data 
and its paramctrization, in the form of a residuals map. We 
then divide, pixel by pixel, the residuals map by the square 
root of the number of counts in the synthetic map, thus 
obtaining a measure of the significance of the disagreement 
in every pixel, which we will call the signal-to-noise ratio, 
S/N, map. 

3. Results 

3.1. Morphology 

To visualize a possible gamma-ray source associated to 
HB21, we produce the S/N map corresponding to the 
null hypothesis, in different energy ranges (Fig. [T]). The 
point spread function^ (PSF) of Fermi /LAT is up to 3° at 
100 MeV energies, im proves to ~ .4° at 1 GeV and achieves 
- 0.1° above 10 GeV (|Randoll2009D . Therefore, at the lowest 
energies (Fig. Ha)) it is not possible to distinguish any struc- 
ture beyond a more or less flat emission extended through- 
out the SNR shell. Above 500 MeV (Fig. [Tb} and above 
1 GeV (Fig. [Tc]) a shell-like structure becomes visible. 

Above 3 GeV (Fig. Ild|) we can identify several struc- 
tures. The most remarkable feature is the NW arc, which 
coincides with the SNR shell and the position of cloud NW. 

3 measured as the 68% event containment radius. 
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(a) S/N map above 100 MeV. Red crosses are the point-like 
sources from 2FGL catalog surrounded by their 95% po- 
sition uncertainty ellipses. Sources highlighted with green 
ellipse are those ones we consider related to HB21 and are 
removed from the model when computing the null hypoth- 
esis map. 




0.2S 0.52 0.78 1 1.3 1.6 1.8 2.1 2.3 2 .6 

(b) S/N map above 500 MeV. The circle used for the spec- 
tral modeling is shown, as well as the divisions discussed in 
tabled] Symbols A, B, C, NW, N, and S mark the position 
of the relevant clouds, mentioned in section [1] The size of 
the markers is not related to the extension of the clouds. 





0.18 0.37 0.56 0.74 . 93 1.1 1.3 1 . 5 1 . 7 l.S 0.1 0.2 0.3 0.4 . 5 0.6 0.7 0.8 0.9 1 

(c) S/N map above 1 GeV. Whit e contours r e prese nt the (d) S/N map above 3 GeV. Green contours depict th e 
large-scale CO distribution from iDame et al.l pOOll ) inte- 4850 MHz radio continuum emission (|Condon et al.lH994l ). 
grated between and -20 km s -1 . Symbols A, B, C, NW, N, 
and S mark the position of the relevant clouds, mentioned 
in section [T] 

Fig. 1: S/N maps above 100 MeV, 500 MeV, 1 GeV, and 3 GeV. The color scale represents signal-to- noise ratio (defined as 
real counts minus model counts divided by square root of model counts) for the null hypothesis. Different multiwavclength 
information is included in each panel. 



There is also a bright spot close to cloud A. This compo- 
nent seems to become more prominent with increasing en- 
ergy. The center of the SNR does not show especially bright 
emission, whereas the S part of the shell presents an en- 



hancement roughly coinciding with cloud S. We show below 
that the emission above 3 GeV is still significant (Tabic [I]). 

To evaluate the morphological properties of the source, 
several scenarios were considered: 
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Tabic 1: TS of the spatial models above different minimum 
energies. N is the number of additional parameters of the 
model with respect to the null hypothesis accounting for the 
spectral index and normalization factor of each additional 
source or component. 



Model 


N 


TSioo 


TS500 


TS1000 


TS3000 


1. 2FGL sources 


10 


959 


610 


245 


48 


2. Circle 


2 


832 


626 


279 


53 


3. 4850 MHz 


2 


780 


613 


275 


42 



1. The four point-like sources from 2FGL coincident with 
HB21; 

2. The 4850MHz map from ()Condon et al.lll994l ). where 
the quasar 3c418.0 was removed from the radio map 
(see discussion below); 

3. A circle of flat emission centered on the catalog position 
of HB21. 

All the templates were rebinned to match the field of view 
and the pixel size of the Fermi /LAT maps. In Table [1] 
we show the TS values for each model in several energy 
regimes, as well as the number of additional parameters in 
each model. The point-like sources model introduces ad- 
ditional degrees of freedom for the flux normalization of 
each of the four sources, plus the spectral indices. 2FGL 
J2043.3+5105 and 2FGL J2046.0+4954 contain an addi- 
tional parameter /3 which allows an energy-dependent spec- 
tral index —a — (3 log (.E/1000GeV). The spectrum of the 
circle and the 4850 MHz templates are described by a 
power-law function at this stage. 

All of the models listed above provide a good description 
of the HB21 field compared to that of the null hypothesis. 
Thus, one or more sources are detected in the HB21 re- 
gion with high significance, even at energies above 3 GeV. 
However, we discard the description by means of point-like 
sources, since it provides the highest TS only at the low- 
est energies (Table [I]), where the broad PSF of Fermi /LAT 
does not allow disentangling any substructure. The model 
with the point-like sources introduces ten degrees of free- 
dom. Therefore we consider that the improvement in TS at 
low energies comes from the greater number of parameters 
and not from the point-like source model being a better 
description of the morphology. With only two degrees of 
freedom, the overall emission can be described well at all 
energies by a circle of radius 1.125° centered on the catalog 
position of the SNR. The radius of the circle is chosen to 
provide the highest TS after varying it from 0.75 to 1.375° 
in steps of 0.125° (one pixel). Initially, we assume a simple 
power-law function as spectral model. 

To further investigate the possibility of a shell-like mor- 
phology we produced the radial profiles of the excess (real 
counts minus model counts, for the null hypothesis) in each 
energy range, and we compared this profile with the one ex- 
pected from models (2) and (3) (Fig.©. In view of the pro- 
files, it is not conclusive whether the gamma-ray emission 
can be related to a shell-like structure. However, given the 
hints obtained so far, we considered this possibility by test- 
ing several ring-shaped models with various inner and outer 
radii. The outer radius is varied within the same range as 
for the optimization of the flat circle, whereas we considered 
inner radii ranging from (circle) to 0.75°, also in steps of 
0.125. The best combination at all energies is 0.125° for the 
inner radius and 1.125° for the outer radius. However, the 



difference in TS is about 0.1 for all energies, which is not 
significant. In addition, we consider that smearing by the 
PSF would prevent such a narrow hole (of only two pixels) 
from being distinguished if used, so we continued to use the 
flat circle for the subsequent analysis. 

We addressed the possibility of having additional point- 
like sources besides the above-mentioned ones. First, the 
NW corner of the radio shell reveals a bright point-like radio 
source due to the presence of the quasar 3c418.0 at z = 
1.6865 (jPaturel et all 12003 ) . We considered the possibility 
that this quasar contributes to the gamma-ray emission. 
To do so, we looked for variability in the GeV signal by 
means of an aperture analysis within a radius of one degree 
around 3c418.0. After several temporal binnings we did not 
find any significant variability in the photon rate around 
the quasar. Therefore, we conclude that if this object has 
any contribution, we cannot disentangle it from that of the 
SNR with the current data set. Secondly, in the 1 GeV map 
(Fig.QxJ) two spots appear SW of the SNR. Comparing with 
the CO large-scale distribution around the SNR we find 
that these spots are roughly coincident with local maxima 
of the gas distribution. Taking this possibility into account, 
we repeated the likelihood analysis adding two sources at 
the position of the gas overdensities. None of these sources 
were significant, whether in the analysis above 100 Me V 
energies or above 1 GeV. 

3.2. Spectral energy distribution 

We divide the considered energy range (100 MeV to 
100 GeV) into twelve bins and compute the spectral en- 
ergy distribution of the whole source by extracting its flux 
in each bin (Fig.©. Only those bins with TS> 10 (~ 3.2cr) 
are shown as spectral points. The last significant bin is the 
one from 3.2 GeV to 5.6 GeV. In addition wc show 95% con- 
fidence level upper limits for the explored energy range, up 
to 100 GeV. The upper limits correspond to the flux pro- 
viding a likelihood value such that 2 A log C = 4. 

Wc notice that the spectrum deviates from a power-law 
function and suggests the presence of a peak at few hundred 
MeV. We test the possibility that the gamma ray emission 
is described by a smoothly broken power law of the form 

Kir) ■ < 2 > 

or a curved power law (log-parabola) 

dN _ ( E \ -"-^STorokv 

~d~E ~ ° V lOOOMcV ) ' ^ ' 

or a power law with a cut-off 

dN »r ( E V 7 E 

= N [ cxp . (4) 

dE VlOOOMeVy E cutoff v ; 

We use the likelihood ratio — 2 \og(C p i / C mo dei) as a mea- 
sure of the improvement of the likelihood fit with respect 
to the simple power law, when different spectral shapes are 
used. We use equations (two extra degrees of freedom), 
[3J and 0] (one extra degree of freedom) , resulting in like- 
lihood ratios of 143, 146, and 130, respectively. Provided 
that the log-parabola introduces an additional parameter 
(/3) to the simple power law, wc conclude that the chance 
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Fig. 2: Radial profile of the excess (real counts minus model counts generated with the null hypothesis model). Red curve 
represents the profile expected from a flat circular emission region of 1.125° radius. Green curve represents the profile 
expected from a gamma-ray emission following the radio continuum emission at 4850 MHz. 



probability of the log-parabola being a better description 
of the spectrum is 1.5 x 10~ 33 . The TS of the flat circle 
(with log-parabolic spectral shape) with respect to the null 
hypothesis is 988, which roughly corresponds to a detection 
significance at the level of 31 standard deviations. The best 
fit parameters are N = (17.5 ± 0.2)KT 12 cm~ 2 s~ 1 MeV~ 1 , 
a = 2.596 ± 0.013, (5 = 0.338 ± 0.008, where uncertainties 
are only statistical. The total energy flux is (1.10 ±0.03) x 
10 -4 MeVcm~ 2 s~ 1 , with a maximum in the spectral energy 
distribution at 413 ± 11 MeV. 

To explore possible spectral differences throughout 
HB21, we divided the circle in three pieces covering 120° 
each and let them acquire different values for the No, ol and 
(3. In this way, segment NE covers clouds A and N; segment 
NW covers cloud NW and segment S covers cloud S. Wc 
used the likelihood ratio — 2 \og{C S egment I ^-circle) ■> where 
^■segment refers to the segmented circle without the tested 
segment. With this likelihood ratio, two thirds of the circle 
may be regarded as part of the null hypothesis, thus provid- 
ing a measure of the significance of the tested segment. All 
three segments contribute significantly to the overall emis- 
sion. The likelihood ratio values of each of these segments 
are summarized in Table [2J along with the flux correspond- 
ing to each segment and their spectral parameters. The 

2 — a 

maximum energy flux is attained at E max = e 2 p GeV. To 



evaluate the uncertainty in E max , the uncertainties in a 
and /3 are taken into account, as well as their covariance. 
Wc find a hint that the segment NW has a softer spectrum, 
peaking at lower energies than the other two regions (see 
also Fig. [3]). 

4. Discussion 

The luminosity between 100 MeV and 5.6 GeV is L = 
(1.34±0.03 sta t) x 10 34 (d/0.8kpc) 2 erg/s. Unless d = 1.7 kpc 
is confirmed, HB 21 belongs to the group of low-luminosity, 
GeV-emitting SNRs, such as Cygnu s Loop (jKatagiri et al.l 
120 111 ) or S147 (jKatsuta et al J 120121 ). which are clearly less 
luminous than the first GeV-e mitting SNRs that w ere dis- 



Abdo et al. 


2010bl) 


Abdo et al. 


2010al) 



, W49B ()Abdo et al.ll2010cD . or Cas A 
have luminosities L > 10 35 erg/s. Also 
the break in energy is found at lower energy in HB 21 than 
in the case of the luminous SNRs. 

Gamma-ray emission from supernova remnants can be 
produced by several non-exclusive mechanisms. Electrons 
and positrons in the supernova remnant shell interact by in- 
verse Compton scattering with ambient photon fields (such 
as infrared starlight) to produce gamma rays. Another 
possibility is that electron-atom or electron-molecule in- 
teractions in a dense medium result in bremsstrahlung 
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Table 2: Spectral analysis of the different circle segments above 100 MeV. 



Segment 


LR 


Flux [Hn B cm -:i s -1 ] (6) 


a (a) 


(3(a) 


E max [GeV] (6) 


Global 





17.5 ±0.2 


2.596 ±0.013 


0.338 ± 0.008 


0.413 ±0.011 


North-East 


90 


4.0 ±0.9 


2.41 ±0.16 


0.33 ± 0.09 


0.54 ±0.14 


North- West 


116 


7.6 ± 1.2 


2.87 ±0.15 


0.32 ±0.07 


0.26 ±0.06 


South 


151 


5.3 ± 1.0 


2.49 ±0.13 


0.39 ±0.09 


0.53 ±0.10 



(a) Spectral parameters of the log-parabolic spectral shape (equation . 
(b) Integral flux and the energy at which the energy flux is maximum. 




Fig. 3: Spectral energy distribution of the gamma-ray emis- 
sion from HB21 modeled as a flat circle of 1.125° radius. 
Red error bar is statistical uncertainty. An additional sys- 
tematic uncertainty of 10% (E < 560 MeV) and of 5% 
(E > 560 MeV) is represented by the black error bar. The 
solid thin black curve is the log-parabola (equation \S§ used 
to model the overall spectrum. Curves with extreme values 
of a and /?, within statistical uncertainty are also shown 
(solid thick black curves). Dotted, dashed, and dash-dotted 
curves are the best spectral descriptions for segments NW, 
S, and NE respectively (see Table [5]). 



radiation. These mechanisms involve electrons acceler- 
ated up to the energies of the observed gamma rays. 
Complementarily, gamma rays can also be produced by 
neutral pion decay, where 7r°s result from the collision 
of accelerated protons (or heavier nuclei) with nucle- 
ons of the ambient gas. The absence of nonthermal X- 
ray emission favors a hadronic origin of the observed 
gamma-ray emission. Moreover, the rapid steepening of 
the spectrum above few GeV is the kind of signature ex- 
pected from the re-acceler a tion of pre-existing cos mic rays 
ijBlandford fc Cowid 119821: lUchivama et~aTI 1201 Oh . where 
high en ergy cosmic rays escape from th e SNR confinement 
region (jZirakashvili fc Ahar onian 20 ftj. 

To check the viability of the lcptonic and hadronic 
scenarios from the energetics point of view, we consider 
the energy from the supernova explosion that is converted 
into accelerated particles, W = L x t. In this expression, 
L is the gamma-ray luminosity and r the characteristic 
cooling time of the dominant accelerated particle type. 
When the gamma-ray luminosity is /ladromc-dominatcd, 
t p ~ 5.3 x 10 7 (7i/lcm~ 3 ) -1 years is the cooling time of 
the accelerated protons as a fun ction of the a mbient pro- 
ton d ensity n (jAharonianl 120041 ) . According to iKoo et all 
(|2001f) . cloud S has density ~ 7000 cm , and central evap- 
orating clouds have densities ~ 4 x 10 4 cm~ 3 . Moreover, 



iTatematsu et alj (|1990l ) quote a density of about 100 cm~ 3 
for cloud A. We computed our own estimate of the average 
density of the region. To estimate the total mass we used 
the CO data from CfA 1.2 m Millimeter- Wave Telescope 
(|Dame et al.lll987t ). We assumed a standard linear relation- 
ship between the velocity integrated CO intensity, IcOj an d 
the molecular hydrogen column density, iV(H2): 

AT(H 2 )/Ico = (1.8 ± 0.3) x 10 20 cm" 2 K- 1 km" 1 s" 1 (5) 

as derived by iDame et al.1 (|200l[ ). This equation yields 
Mco/Mq = 1 200Sco d-lpc where dkpc is the distance 
to the cloud in kpc, and Sco the CO emission inte- 
grated over velocity and the angular extent of the cloud 
in K km s _1 arcdegree 2 . We conclude that there are as 
many as 12000 solar masses of molecular gas coinciding 
with the gamma-ray emission in the velocity range from 
-20 to Okm/fl Assuming that this gas is in a spherical 
volume of 25 pc diameter we obtain an average density 
of about 60 protons per cubic centimeter, and therefore 
W p ~ 4 x 10 47 crg. Instead, in case the gamma-ray lu- 
minosity is Zepfomc-dominated, the cooling time due to 
brcmsstrahlung interactions depends on the ambient den - 
sity as Tb rems ~ 4 x 10 7 (n/cm 3 ) -1 ([Aharonia n| 12001 . 
which leads to a live time T^ rems ~ 7 x 10 5 years and re- 
quires an energy budget similar to the hadronic interac- 
tions, W e ~ 3 x 10 47 erg. Alternatively, the cooling time 
of electrons due to synchrotron loss es is T sync ~ 1.3 x 
10 10 ( J B/l^G)- 2 (£: e -/lGcV)- 1 years dGaisser et al.lll998lh 
Therefore, an E e - = 3 GeV electron in a typical inter- 
stellar magnetic field of 3 /j,G would have a live time of 
T S ync ~ 5 x 10 8 years. Only with a magnetic field as 
high as 80 /iG, could the synchrotron losses become effi- 
cient enough to make them comparable to the losses due 
to bremsstrahlung interactions. Although detailed model- 
ing is needed to confirm these considerations, our order-of- 
magnitude estimates favor the view that the high-density 
regions surrounding the SNR determine the dominant emis- 
sion process. This is supported further by the coincidence of 
the gamma-ray bright regions with dense clumps of molec- 
ular gas (Fig. [I}. In any case, energy budget considerations 
are not sufficient to distinguish between electron-dominated 
and proton-dominated scenarios under the assumption that 
the supernova explosion has a typical energy release of 
about 10 51 erg. 

Northeast of the circle there is a spot that becomes 
bright with increasing energy (Fig.QJ. We discuss the possi- 
bility that this comes from a hard gamma-ray emission com- 
ponent related to cloud A. This cloud most likely already 

4 iBvun et al.1 lj2006T ) show that the majority of the emission 
is concentrated in the range -20 to km/s. Although there is 
some emission at higher velocities, at about 3 km s _1 , telluric 
CO emission corrupts the spectra. To avoid this issue, and for 
simplicity, we took the velocity range -20 to km/s. 
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existed at the time of the SNR explosion, it is not shocked, 
and its v elocity is not significan tly different fr om that of 
the w all (|Tatematsu et al.lll990D . According to iKoo et al.l 
(|200lF) . this cloud does not present broad line emission, 
so there is no direct evidence of interaction with HB210. 
Moreover, cloud A coincides with a concavity of the SNR 
shell. The coincidence of cloud A with this feature in the 
radio continuum emission suggests that the overtaking of 
this cloud leads to retardation of the shock front in com- 
parison to the surrou n ding although high-resolution CO 
maps from IKoo et al~l (|2001h did not find such evidence. 
There is still the possibility that the shock is dissociative, 
and molecules have not been reformed (thus being omit- 
ted in the search of broad line emission regions), but the 
shock velocity of 20 km s _1 observed by Koo is in princi- 
ple not enough for molecule dissociation (which typically 
requires 25-50 km s _1 ). Finally, the infrared emission of 
cloud A is caused by dust heated by the ambient radiation 
fields in their surface, whereas in clou ds N and S it is due 
to li nes from sh o ck-ex cited molecules (iKoo et al.ll2001l) . In 
fact, IKoo et al.l (|200lD mentions a diffuse component con- 
necting clouds N and A. All in all, there is evidence that 
cloud A is different to other molecular clouds in the vicinity 
of HB21. As mentioned in section [1] there is also the sug- 
gestion that cloud A could be in the foreground of HB 21, 
related to the Cyg OB7 association (i.e. at 0.8 kpc), whereas 
HB21 could be in the background, at 1.7 kpc or more. If 
we assume that the gamma-ray brightness of cloud A is 
due to runaway protons from HB21, we can estimate the 
maximum distance between the two objects by the relation 
R d = VWi, where D ~ lO^cn^s" 1 (jGabici et al.ll2009l 
equation 11) is the diffusion coefficient of cosmic rays pro- 
tons of 10 GeV (which originate 1 GeV gamma rays). In this 
case the separation between cloud A and HB21 would be 
roughly 50 pc. We note that D is completely unknown, and 
this distance is not to be taken as a measure of the separa- 
tion between the cloud and the SNR, but only a suggestion 
that cloud A is indeed close to the SNR. We cannot con- 
clude the same about the two other clouds ( B and C) on 
the eastern rim detected bv lTatematsu et all (|1990t i. These 
are dark in gamma rays and th erefore they may b e well in 
the foreground, as suggested by iBvun et al.l (|2006l ). 

Finally, we consider it interesting that the part of the 
circle related to cloud NW presents a spectral break at lower 
energies than the rest of the SNR. From the E > 3 GeV map 
(Fig. Ild[) . we see that the majority of the emission in this 
region comes from a bow-shaped structure that resembles 
the SNR shell itself. We suggest that the gamma-ray emis- 
sion from this region originates in the shell itself, or from a 
molecular cloud (probably cloud NW) that was overtaken 
at an earlier stage than those producing the emission in the 
remaining two thirds of the remnant. 

5. Conclusions 

The analysis of 3.5 years of public Fermi /LAT data leads to 
a clear detection of an extended source of gamma rays coin- 
cident with the SNR HB 21. Our own method for visualizing 
the morphology allows us to describe the source as a flat 
circle with uniform emission, although we resolve a clumpy 

5 The lack of broad line emission is also the case for 
3c391, where the interaction with molecular clouds is certain 
(jReach et al.ll2002l ). 



structure at the highest energies. The spectral analysis re- 
veals a peak of the gamma-ray emission at 413 ±11 MeV. 
However, we conduct a dedicated spectral analysis for three 
regions, which show particularly bright emission above 
500 MeV coincident with dense molecular hydrogen clumps. 
We find indications that the gamma-ray emission peaks at 
somewhat lower energies (0.26 ± 0.06 GeV) in the NW re- 
gion of the shell (coincident with cloud NW and the shell 
itself), whereas in the NE (coincident with clouds N and 
A) and the S regions (coincident with cloud S) the peak is 
found around ~ 0.5 GeV. 

Such spectral breaks are expected in middle-aged SNRs 
like HB 21 due to the escape of CRs from the confinement 
region. In HB 21 the break occurs at lower energies than in 
other similar objects. This could be related to the proximity 
of very dense molecular clouds to the supernova explosion 
progenitor, which c ould have slowed down the shock rapidly 
(jOhira et al.ll201lF) . The spectrum from the region related 
to the bow-shaped emission in the NW peaks at lower en- 
ergies than the spectrum from clouds N and S. This fact 
matches the understanding that the most energetic par- 
ticles related to the SNR may have already escaped the 
NW region, due to a smaller confinement volume, or be- 
cause cloud NW was shocked at an earlier stage than clouds 
N or S. Moreover, in the S/N map above 3GcV, we see 
how a spot coincident with cloud A that becomes relatively 
brighter than the other clouds with increasing energy. We 
suggest that this is because cloud A is separated from the 
SNR, and it is now reproducing the spectrum of the SNR 
at an earli er stage due to particles tha t diffused away from 
the shock (jAharonian fc At ovanll 19961 ). 

Therefore, given that HB 21 appears as an extended 
object even for gamma-ray telescopes, it provides the op- 
portunity to observe the production and diffusion of accel- 
erated particles (most likely protons), from the SNR shell 
to distant molecular clouds acting as targets. 
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